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We study excitations in a normal Fermi liquid with a local scalar interaction. Spectrum of bosonic 
scalar-mode excitations is investigated for various values and momentum dependence of the scalar 
Landau parameter /o in the particle-hole channel. For /o > 0 the conditions are found when the 
phase velocity on the spectrum of the zero sound acquires a minimum at a non-zero momentum. For 
— 1 < /o <0 there are only damped excitations, and for fo < — 1 the spectrum becomes unstable 
against a growth of scalar-mode excitations (a Pomeranchuk instability). An effective Lagrangian for 
the scalar excitation modes is derived after performing a bosonization procedure. We demonstrate 
that the Pomeranchuk instability may be tamed by the formation of a static Bose condensate of 
the scalar modes. The condensation may occur in a homogeneous or inhomogeneous state relying 
on the momentum dependence of the scalar Landau parameter. Then we consider a possibility of 
the condensation of the zero-sound-like excitations in a state with a non-zero momentum in Fermi 
liquids moving with overcritical velocities, provided an appropriate momentum dependence of the 
Landau parameter fo(k) > 0. 
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I. INTRODUCTION 

The theory of normal Fermi liquids was built up by 
Landau pj, see in textbooks M- The Fermi liquid ap¬ 
proach to the description of nuclear systems was devel¬ 
oped by Migdal fo. (6|. In the Fermi liquid theory the low- 
lying excitations are described by several phenomeno¬ 
logical Landau parameters. Pomeranchuk has shown in 
Ref. (zl that Fermi liquids are stable only if some inequal¬ 
ities on the values of the Landau parameters are fulfilled. 

In this work we study low-lying scalar excitation modes 
(density-density fluctuations) in the cold normal Fermi 
liquids for various values and momentum behavior of the 
scalar Landau parameter fo in the particle-hole channel. 
We assume that an interaction in the particle-particle 
channel is repulsive and the system is, therefore, stable 
against pairing in an s-wave state. An induced p-wave 
pairing to be possible at very low temperatures T < T c>p , 
see Refs. 5, can be precluded by the assumption that 
the temperature of the system is small but above T c , p . 

For /o > 0 the conditions will be found when the 
phase velocity of the spectrum possesses a minimum at 
a non-zero momentum. This means that the spectrum 
satisfies the Landau necessary condition for superfluid¬ 
ity. As a consequence this may lead to a condensation 
of zero-sound-like excitations with a non-zero momentum 
in moving Fermi liquids with the velocity above the Lan¬ 
dau critical velocity [9j. Similar phenomena may occur 
in moving He-II, cold atomic gases, and other moving 
media, like rotating neutron stars, cf. Refs. mm- ^ 
— 1 < fo < 0 excitations are damped, for f 0 < — 1 the 
spectrum is unstable against the growth of zero-sound¬ 
like modes and hydrodynamic modes. Up to now it was 
thought that for /o < — 1 the mechanical stability condi¬ 
tion is violated that results in exponential buildup of the 


density fluctuations. In hydrodynamic terms the condi¬ 
tion /o < — 1 implies that the speed of the first sound 
becomes imaginary. This would lead to an exponential 
growth of the aerosol-like mixture of droplets and bub¬ 
bles (spinodal instability). In a one-component Fermi liq¬ 
uid spinodal instability results in a mixed liquid-gas like 
stationary state determined by the Maxwell construction 
if the pressure has a van der Waals form as a function 
of the volume. In the isospin symmetric nuclear matter 
(if the Coulomb interaction is neglected) the liquid-gas 
phase transition might occur [161I17| for the baryon den¬ 
sities 0.377,0 < n £ 0.7no, where no is the nuclear satura¬ 
tion density. In a many-component system a mechanical 
instability is accompanied by a chemical instability, see 
Ref. 18|. The inclusion of the Coulomb interaction, see 
Refs, jl9|,[20|, Fads to a possibility of the pasta phase in 
the neutron star crusts for densities 0.3no ;$ n 0.7no. 

The key point of this work is that we suggest an al¬ 
ternative description of unstable zero-sound-like modes 
which might be realized at certain conditions. We 
demonstrate that for fo < — 1 instability may result in 
an accumulation of a static Bose condensate of the scalar 
field. The condensate amplitude is stabilized by the re¬ 
pulsive self-interaction. The condensation may occur in 
the homogeneous either in inhomogeneous state depend¬ 
ing on the momentum dependence of the Landau param¬ 
eter fo(k). In the presence of the condensate the Fermi 
liquid proves to be stable. 

The work is organized as follows. In Sect. [II] we study 
spectrum of excitations in a one-component Fermi liquid 
in the scalar channel in dependence on fo(k). In Sect. Mil 
we bosonize the local interaction and suggest an effec¬ 
tive Lagrangian for the self-interacting scalar modes. In 
Sect. [TV] we study Pomeranchuk instability for / 0 < — 1 
and suggest a novel possibility of the occurrence of the 
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static Bose-condensate which leads to a stabilization of 
the system. In Sect.[V]we consider condensation of scalar 
excitations in moving Fermi liquids with repulsive inter¬ 
actions. Concluding remarks are formulated in Sect. IVII 


II. EXCITATIONS IN A FERMI LIQUID 
A. Landau particle-hole amplitude. 


Consider the simplest case of a one-component Fermi 
liquid of non-relativistic fermions. As discussed in the 
Introduction, we assume that the system is stable against 
pairing. The particle-hole scattering amplitude on the 
Fermi surface obeys the equation @-0, [§] 

T p h(n', n; q) = Y u (n' ,n) 

+ {ft',n”) C ph {n"; q) f ph (n",n; q)) H " , (1) 

where n and n ’ are directions of fermion momenta before 
and after scattering and q = (w, k) is the momentum 
transferred in the particle-hole channel. The brackets 
stand for averaging over the momentum direction n 

<->" = /tt(->. < 2 > 

and the particle-hole propagator is 

+oo +oo 

£ P h(n; q) = J ^~J g (pf+) G(pf-) , (3) 

— oo 0 


where we denoted pf± = (e ± uj/ 2 ,pFn ± k/ 2 ) and pp 
stands for the Fermi momentum. The quasiparticle con¬ 
tribution to the full Green’s function is given by 


G(e,p) = 


e - + i 0 signe 


(4) 

^ 2m* F • [ ’ 


Here m F is the effective fermion mass, and the param¬ 
eter a determines a quasiparticle weight in the fermion 
spectral density, 0 < a < 1, which is expressed through 
the retarded fermion self-energy S F (e,p) as a” 1 = 1 — 
(95REp/9e)o. The full Green’s function contains also a 
regular background part G reg , which is encoded in the 
renormalized particle-hole interaction T“ in Eq. 0- 
The interaction in the particle-hole channel can be 
written as 


f“(n', n) = Tq (n'n) CTqCTo + T" (n'n) (a'a). (5) 

The matrices Uj with j = 0,..., 3 act on incoming 
fermions while the matrices er' act on outgoing fermions; 
do is the unity matrix and other Pauli matrices 0 - 1 , 2,3 are 
normalized as Trcqerj = 2(5,j. We neglect here the spin- 
orbit interaction, which is suppressed for small trans¬ 
ferred momenta q <C pf- The scalar and spin amplitudes 


in Eq. © can be expressed in terms of dimensionless 
scalar and spin Landau parameters 


f(n',n) = a 2 N 0 Y%(n',n ), 

g(n',n) = a 2 N 0 T“{n',n) , (6) 

where the normalization constant is chosen as in applica¬ 
tions to atomic nuclei [b|, HH [HJ with the density of states 
at the Fermi surface, Nq = N(n = no), taken at the nu¬ 
clear saturation density n 0 and N (n) = m F( n )P F ( n ) _ Such 
a normalization is at variance with that used, e.g., in 
Refs. HH3- Their parameters are related to ours defined 
in Eq. © as / = N f /N 0 and g = N g/N 0 . 

The Landau parameters can be expanded in terms of 
the Legendre polynomials Pi(n ■ n r ), 

f{n',n) =^2fiPi(n-n') , (7) 

1 

and the similar expression exists for the parameter g. 
The Landau parameters /o,i, go,i can be directly related 
to observables [1]. For instance, the effective quasiparti¬ 
cle mass is given by m 

—= 1 + a 2 NT 0 (cos 9) cos 0 = 1 + , (8) 

ttif 3 

where the bar denotes the averaging over the azimuthal 
and polar angles. The positiveness of the effective mass 
is assured by fulfillment of the Pomeranchuk condition 
fi > —3. Note that the traditional normalization of the 
Landau parameters (O depends explicitly on the effective 
mass to* through the density of states N. Therefore, it 
is instructive to rewrite Eq. ([5]) using the definition in 
Eq. © 


TOp _ _1_ 

m F “ 1 - \-niFj ■ 
3 (no) J 


( 9 ) 


From this relation we obtain the constraint /1 < 

3TOp(no)/TOF for the effective mass to remain positive 
and finite; otherwise the effective mass tends to infinity 
in the point where /1 = 3 m F (no)/ tuf■ Thus, for the sys¬ 
tems where one expects a strong increase of the effective 
mass, the normalization 0 of the Landau parameters 
would be preferable. A large increase of the effective 
fermion mass may be a s ign of a quantum phase transi¬ 
tion dubbed in Refs. [23l . [+it as a fermion condensation. 
The latter is connected with the appearance of multi- 
connected Fermi surfaces. As demonstrated in Ref. j25j if 
such a phenomenon occurs in neutron star interiors, e.g., 
at densities close to the critical density of a pion con¬ 
densation, this would trigger efficient direct Urea cooling 
processes. On the other hand, right from the dispersion 
relation follows that 

to * _ 1 - (cMZ*/de) F 

tof 1 + (95REp/i9e°)F ’ 
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where e° = p 2 /2mp. Thereby, the Landau parameter jy 
can be expressed via the energy-momentum derivatives 
of the fermion self-energy. 

Below we focus our study on effects associated with the 
zero harmonic fo in the expansion of Tq as a function of 
(nn'). Then the solution of Eq. © is 

T ph (n',n;q) = T ph ,o(g) 4 a 0 + T ph ,i(g) (ct'ct) , 
Tph ’° (1)(9) = l/T“ { 1 ) -(C ph (n-, q )) H - (U) 


The averaged particle-hole propagator is expressed 
through the Lindliard’s function 

(£ ph (n; q))n = - a 2 N<f >(—) , (12) 

v vf k pf 


where 


<&(s, x) 


Z--1 Z- + 1 

—7-- in- 

4 (z+ - z-) z- - 1 

4-1 ,f t +i + i 

4 (z+ — Z-) z+ — 1 2 


(13) 


Here and below we use the dimensionless variables z± = 
s ± x/2 , s = ui/kvF , a; = fc/pp. For real s the Lindhard’s 
function acquires an imaginary part 


3$(s,a;) = 
For s>l, 


f s , 0<s<1-f 
iy(l-^ 2 ) , 1-f <s< 1 + f ■ (14) 
0 , otherwise 


$(s,4 « -1/(3 2+*_). (15) 

For x -C 1 the function $ can be expanded as 

$(s, x) ps 1 + ^ log -i-—- 2 • (16) 

V ; 2 S s +1 12 (s 2 — l) 2 ^ 

and if we expand it further for s<Cl we get 


<f>(s, x) 


7r 9 x s 

1 + — s 2 — — —-— 


12 


-So- < 171 


At finite temperatures the function $ should be re¬ 
placed by the temperature dependent Lindhard function 
calculated in Ref. [26|. Generalization of expan¬ 
sion m for low temperature case (T <C cf) is given 
by 


/ v 

$ T (s,X,t) = $(s,x) 4 - — t 2 J , 


(18) 


where t = T/e p and ep = p|/2?Up is the Fermi energy. 
For high temperatures (T e p) we have 


$T(0,X,t) 


—(i~ — —2=). 
3 1\ 6 1 3 V2tt t 3 ' 


(19) 


The amplitudes (fill) possess simple poles and logarith¬ 
mic cuts. For the amplitude Tp^o the pole is determined 
by the equation 


1 

To 


-$( 


lo k 
vf k ’ pf 


( 20 ) 


Similar equation exists in ^-channel (with replacement 
/o “^ 9o- Analytical properties of the solution (EOl) have 
been studied in [2?| . 

Expanding the retarded particle-hole amplitude 
T* 0 (q) near the spectrum branch 


2 u(k)V 2 (k) 


Tph ’ o(g ) ^ (co + i0) 2 -uj 2 (k) ’ 


V~ 2 (k) = a 2 N 


d<Z> 

dJ 2 


uj(k) 


( 21 ) 


we identify the quantity D r (lo, k) = [(w + *0) 2 — w 2 (fc)] _1 
as the retarded propagator of a boson with the dispersion 
relation ui = ui(k) and the quantity V(k) as the effective 
vertex of the fermion-boson interaction. 

For the neutron matter the parameters / = f nn , 
g = gnn are the neutron-neutron Landau scalar and 
spin parameters. Generalization to the two component 
system, e.g., to the nuclear matter of arbitrary isotopic 
composition is formally simple [§]. Then the amplitude 
should be provided with four indices ( nn , pp, np and 
pn). However, equations for the partial amplitudes do 
not decouple. For the isospin-symmetric nuclear sys¬ 
tem with the omitted Coulomb interaction the situation 
is simplified since then f nn = f pp and f np = f pn . In 
this case one usually presents Tp in Eqs. © and © as 
r£ = (/ + /W)/« 2 fVo, and ry = (g + g'f'r)/a 2 N 0 , 
where r,; are isospin Pauli matrices. Quantities / and f 
(and similarly g and <?') are expressed through f nn and 
fnp as f 2 (fnn T fnp ) and f — 2 (fnn fnp)- The am¬ 
plitudes of these four channels (/, /', g , g') decouple, cf. 
Ref. & The excitation modes are determined by four 
equations similar to Eq. (l2l)l) . now with 2/o, 2/g, 2go 
and 2 g' Q . Presence of the small Coulomb potential mod¬ 
ifies the low-lying modes determined by Eq. (|20l) (now 
with jo, jo, go and g' 0 ) only at low momenta (28|. We 
will ignore the Coulomb effects in our exploratory study. 

We discuss now the properties of the bosonic modes 
for different values of the Landau parameter /q. 


B. Repulsive interaction (/o > 0). Zero sound. 

For /o > 0 there exists a real solution of Eq. m such 
that for k —> 0 the ratio u> s (k)/k tends to a constant. 
Such a solution is called a zero sound. The zero sound 
exists as a quasi-particle mode in the high frequency limit 
ui l/r co i, where r co i oc £f/T 2 is the fermion collision 
time. In the opposite limit it turns into a hydrodynamic 
(first) sound [29(. At some value /ci im < pp the spectrum 
branch enters the region with > 0, and the zero sound 
becomes damped diffusion mode. 
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We search the solution of Eq. (l20l) in the form 

u>s(k) = k vf s(x), ( 22 ) 

where s(x) is a function of x = k/pF which we take in 
the form s(x) ~ sq + s 2 x 2 + S 4 x 4 . The odd powers of x 
are absent since is an even function of x. 

Above we assumed that / depends only on fin'. The 
results are however also valid if the Landau parameter / 
is a very smooth function of x 2 . From now we suppose 

fo(%) ~ foo + fo 2 X 2 , (23) 


where the parameter /02 is determined by the effective 
range of the fermion-fermion scattering amplitude and 
expansion is valid provided |/oo| | fm I for relevant 

values x < pf- According to Ref. [2l|, in the case 
of atomic nuclei (n ~ n 0 ), /02 = ~ f00^1^1?/%, and 
0.5 < r e ff ~ lfm, as follows from the comparison with 
the Skyrme parametrization of the nucleon-nucleon in¬ 
teraction and with the experimental data. For isospin- 
symmetric nuclear matter /oo > 0 for n ^ no, foo < 0 for 
lower densities, and in a certain density interval below 
n o, fo < — 1, cf. Refs. mil- In the purely neutron 
matter one has — 1 < foo <0 for n f: no |31| . 

The constant term, so, follows from Eq. © 


1 + foo _ £0 j Sp + 1 
foo 2 So -1 

For foo -C 1 the solution of Eq. (l24l) is 

so = 1 + C[1 + (4 + 5/oo)(C'/2/oo) 

+ (24 + 52/qo + 29/o 0 )(C/2/oo) 2 ] , 


(24) 


(25) 


where C = (2/e 2 ) exp(— 2// 0 o). For f 0 0 < 1.8 this ex¬ 
pression reproduces the numerical solution with a percent 
precision. In the oppos ite limit f 00 1 the asymptotic 

solution is so = \/ /oo/3. The coefficients S 2 and S 4 follow 
as 



f02 

1 d 2 fi> 


|"<9<I> 

n 

S2 = 

k 2 o 

2 dx 2 

so,0_ 

V ds 

s 0 ,0- 


(26) 


So(o + / 02 )(s§ - 1) 
foo (1 + foo — So) 


= /oo/[f2(Sg - l) 2 ] . 


S4 — So/oo- 


2 1 l+5 s n — S0 (-^q — l)spS 2 _ fo2 ( „2 
S 2 "l 240 (sq —l) 2 ft\ S 0 


(sq 1)(1 T foo Sq) 



The numerical solutions of Eq. (BUll are shown in Fig. [T] 
by solid lines for various values of the Landau parameters 
foo and /o 2 - The zero-sound solution exists for s 0 > 0. 
We check that the inequality 1 + foo — Sq >0 holds for 
foo > 0. The quadratic approximation for the spectrum 
u) ~ uf k [so T S 2 X 2 ] is demonstrated by dashed lines. 
Only the quasi-particle part of the spectrum is shown 
with ui s (k)/ (vf k) > 1 + k/(2pF). This inequality holds 
for k/pF < fcmax/PF = (1 - \/1 - 16s 2 (s 0 - l))/(4s 2 ) if 



FIG. 1: Spectrum of the sound modes defined by Eq. (© for 
various values of the Landau parameters /oo and /o 2 - Solid 
lines show the results of the numerical solution, dashed lines 
are the quadratic expansions with the parameters so and S 2 
given by Eqs. (1241) and (1261) . The values of foo are shown in 
labels. Lines with (without) stars are calculated for /02 = 0 
(/02 = —1.5a), where a is defined in Eq. (l26l) . The dash- 
dotted line shows the border of the imaginary part of the 
Lindhard’s function m- Below this line the frequency is 
complex. 


S 2 < l/[16(so — !)]■ For larger k the Lindhard func¬ 
tion becomes complex for continuation of the branch 
w = uj s (k) and the solution of Eq. (l20ll acquires an 
imaginary part. For S 2 > l/[16(so — 1)], the function 
oj s [k) does not enter in the region of the complex Lind¬ 
hard function. However, for positive S 2 there is another 
source of the mode dissipation related to the decay of 
one mode’s quantum in two quanta. The latter pro¬ 
cess is allowed if the energy-momentum relation lo s ( p ) = 
u) s (p') +uj s {\p—p'\) holds, that is equivalent to the rela¬ 
tion 1 — cos (f> = 3s2/(soPf) ) where cos= {pp')/(pp'), 
which is fulfilled if S 2 > 0. 

In Fig. |T] we demonstrate first the case fo = const, i.e., 
/02 = 0 . As we see, in this case S2 > 0 for any f 0 0 > 
0 . We see that the quadratic approximation coincides 
very well with the full solution. Then we study how the 
spectrum changes if the parameter /02 is taken nonzero. 
As we conclude from Eq. © the coefficient S2 can be 
negative for /02 < —a. The latter implies: /02 ^ —10 
for foo = 1; fo2 < -2.6 for f 00 = 2; and f 02 Z -1.7 
for foo = 3. In Fig. |T] by solid curves and dashed curves 
marked with stars we depict the zero-sound spectrum 
for fo2 < —a, being computed following Eq. ® and 
within the quadratic approximation, respectively. We 
chose here /02 = — 1 . 5 a. The quadratic approximation 
for s(x) works now worse for momenta close to fc max and 
the next term S4X 4 should be included. We note that 
the parameter S4 is positive in this case and the function 
s(x) has a minimum at x m i n = \f\s 2 \/{ 2 sf). In the point 
k = ko corresponding to the minimum of io s {k)/k the 
group velocity of the excitation u gr = dui s /dk coincides 
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FIG. 2: Coupling constant of zero-sound modes and fermions 
as a function of the momentum k for various values of the 
Landau parameters /oo and / 02 , calculated from Eq. (1211) . All 
curves end at the limiting value of a; max = fcmax/pF at which 
the spectral branch enters into the region of complexity of the 
<1? function. 



k/ Pp 

FIG. 3: Solid lines show the purely imaginary solutions of 
Eq. (1311) for various values of the Landau parameter /oo < 
0. Dashed line indicates positions of the maxima of spectra 
u >d for /oo < —1. Dash-dotted lines show the approximated 
spectra tud(fe) ~ wf k(sdo+Sd 2 x 2 +Sd 4 X 4 ) with the parameters 
from Eqs. (I5ll) and (1551) . 


with the phase one r) p h = io s /k. The quantity u> s {ko)/ko 
coincides with the value of the Landau critical velocity 
for the production of Bose excitations in the superfluid 
moving with the velocity u > u^. 

The effective vertex of the boson-fermion coupling HID 
is shown in Fig. [2] as a function of x for several val¬ 
ues of /oo and /o 2 - We plot it in the range 0 < x < 
Xmax = fcmax/pF, where the quasiparticle branch of the 
zero-sound spectrum is defined. For /02 = 0 the vertex 
V(k) is shown by solid lines, for /02 < 0 by dash-dotted 
lines. We see that V 2 (a; max ) is smaller for / 02 < 0 than 
for /02 = 0. For small x one can use the analytical ex¬ 
pression 


a 2 N 


V 2 (x) 

Vf PF 
S2Xf? 


r<9$ 


L ds 

s(x),x- 


00 


1 + 


4as 2 a; 2 


ot + S 02 L 0 . + /02 


( 6 s 2 



(28) 


which works well for k < 0.75fc max . For k —> 0 we get 
a 2 NV 2 {k ) ~ s 0 (s§ ^ f) 2 /oo^F/(l + /oo - «o)- 


More generally, for purely imaginary s = i s, s £ R, 
the Lindhard function can be rewritten as 


<I>(zs, x) = $(s, x) = — \s 2 - jX 2 + ll log 
Ax * J 


s s 

arctan-,-p arctan- 


1 - 2 X 


1 + 


~s 2 + (\x+l) 2 
s 2 + {\x-l) 2 
1 — 7rs 

2 ' 

(30) 


We note that the function arctan(s/[l — \ x]) should be 
continuously extended for x > 2 as arctan(s/[l — \x}) + 
7 rsigns . Solutions of equation 

l/l/ool = $(*s,x) (31) 

for 0 > /oo > — 1 are depicted in Fig. [3] for / Q0 = 
—0.2; —0.5; —0.9. These solutions are damped (—iw < 0). 
For /02 ^ 0 (at the condition 0 > fo(k) > —1), the 
damped character of the solutions does not change. 


C. Moderate attraction, —1 < /oo < 0. Diffusons. 


D. Strong attraction, foo < — 1. Pomeranchuk 
instability. 


Assume /02 = 0. For —1 < f 00 < 0, Eq. (l20l) has only 
damped solutions with Jis < 1 and 3s < 0. In the limit 
of s,x 1, using the expansion 113 we easily find the 
analytic solution 


Ud(k) 

kv f 


*s d (a;) 


.(2 1 — | /oo | 

1 v V l/oo | 


x 2 

12 


(29) 


which is valid for 1 — |/oo| 1 and x <C 1. We see 

that the solution we found is purely imaginary and its 
dependence on x is very weak. 


Now we consider the case of a strong attraction in the 
scalar channel /oo < — 1. We continue to assume /02 = 0 
for certainty. The Pomeranchuk instability for foo < — 1 
manifests itself in the negative compressibility which is 
found from the relation for the variation of the fermion 
contribution to the chemical potential [ 2 ] 


K = 


pK 71 0 ) 

3 TOp(no) 


(1 + /oo), 


(32) 


being negative for foo < — 1. 
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Solutions of Eq. (1311) for f 00 < — f are shown in Fig. [3] 
We see that there is an interval of k, 0 < k < kd, where 
—iojd > 0 and, hence, the mode is exponentially growing 
with time. This corresponds to the Pomeranchuk insta¬ 
bility in a Fermi system with strong scalar attraction, 
see Ref. [I]. For k > kd we get —*w d (fc > fed) < 0, and 
the mode becomes stable again. The momentum kd is 
determined from Eq. (1311) for s = 0: 


1 __ x\ - 4 x d + 2 
foo 8z d ° S x d - 2 


1 

2 ’ 


(33) 


Xd = kd/pF- The instability increment — iw d has a max¬ 
imum at some momentum k m . The locations of maxima 
for different foo values are connected by dashed line in 

Fig. El 

The spectrum of the unstable mode can be written for 
x « 1 as Sd(a:) ~ Sdo + Sd 2 % 2 + s d 4 x i , where the leading 
term is determined by the equation 

7T Zf 

- - arctans d0 = — , z f = 1 - 1/|/ 00 | • (34) 

z SdO 

The subleading terms are equal to 


Sd2 

Sd4 


s d 0 


12 ( 1 + s do)( z /( 1 + s do)^ s do) 


1 - 5f|p 

20(1 + a 2 d0 ) 


;Sd2 — 


4s d oS 2 


1 + s d0 


d2 - 12s 2 


d2 ■ 


(35) 


This approximation of the spectrum is illustrated in 
Fig. El by dash-dotted lines. The approximation works 
very well for —1.5 < foo < —1 (dash-dotted lines and 
solid lines coincide) but becomes worse for smaller foo- 
For a slightly subcritical case, where 0 < Zf <C 1, we 
obtain 


uid(k) s=s i-v F k [Zf - fc 2 /(12p|)]. (36) 

7 r 

The function —ico d has a maximum at k m = 2pFyfzJ 
equal to 


- iui m = (8/ 3tt)vfPf z 3 / 2 . (37) 


E. Spinodal instability. 


limit to the collision-less regime of the zero sound, i.e. 

For the van der Waals equation of state the compress¬ 
ibility and the square of the first sound velocity prove 
to be negative in the spinodal region. Thereby excita¬ 
tion spectrum is unstable and in simplest case of ideal 
hydrodynamics the growing mode is as follows 32 1 


—iuj = k\J\u 2 — ck 2 1 


(39) 


where c is a coefficient associated with the surface ten¬ 
sion of the droplets of one phase in the other one. Note 
that the viscosity and thermal conductivity may delay 
formation of the hydrodynamic modes. Maximum in k 
yields 


-iuim = v F PF(\foo\ - l)/(6?7i FV / c) • (40) 


The rate of the growing of the spinodal mode decreases 
with increase of the surface tension of the droplets, 
whereas the growing rate of the collision-less mode does 
not depend on the surface tension. For 


^ 2 |/oo | 3/2 
256m|(|/ 00 | - 1) 


(41) 


the zero-sound-like excitations (1361) would grow more 
rapidly than excitations of the ideal hydrodynamic mode 

m- 

For isospin symmetric nuclear matter the spinodal re¬ 
gion in the dependence P(l/n) exists at nucleon densi¬ 
ties below the saturation nuclear density, see Ref. 53- 
In multi component systems with charged constituents, 
like neutron stars, the resulting stationary state is the 
mixed pasta state, where finite size effects (surface ten¬ 
sion and the charge screening) are very important, cf. 
Refs. 0, [33], contrary to the case of the one component 
system, where the stationary state is determined by the 
Maxwell construction, cf. Ref. E 3 - The liquid-gas phase 
transition may occur in heavy-ion collisions. A nuclear 
fireball prepared in a course of collision has a rather small 
size, typically less or of the order of the Debye screening 
length. Therefore, the pasta phase is not formed, as in a 
one-component system. 


For foo < —1 not only the compressibility is negative 
but also the square of the first sound velocity j2J 


III. BOSONIZATION OF THE LOCAL 
INTERACTION 


u 2 = 


Pi 


dP 


dp 3m F Jnp 


-(l + /oo)<0. 


(38) 


Here P denotes the pressure and p is the mass density. 
Note that in the limit T —> 0 the isothermal and adiabatic 
compressibilities coincide, whereas for T ^ 0 the differ¬ 
ence becomes substantial, see Ref. [32]. Also note that 
the first sound exists in the hydrodynamical (collisional) 
regime, i.e. for w <C oc T 2 , which is the opposite 


The description of a fermionic system with a con¬ 
tact interaction can be equivalently described in terms of 
bosonic fields <p q = here we denote q = (w, k) 

and p = (e,p). In terms of the functional path inte¬ 
gral the transition to the collective bosonic fields can 
be performed with the help of a formal change of vari¬ 
ables by means of a Hubbard-Stratonovich transforma¬ 
tion EJ, |35|. After this transformation the effective Eu¬ 
clidean action for the system with a repulsive interaction, 
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here fo > 0, can be written within the Matzubara tech¬ 
nique in terms of the bosonic fields as 


Then making the replacement (1461) in the fermion scat¬ 
tering amplitude CD we obtain 


S [4>] = Y 

9 



Tr log 1 


Gi(j) q ■ 
0-1’ 


— (a 2 N) Vow 


( 42 ) T ^' k) = « 0 )-W)-£gVV) 55 . 


(47) 


where G and (j) are infinite matrices in fre¬ 
quency/momentum space with matrix elements 

= ^p,p'G(c : p^ and [^]g,g / = ^ 9 , 9 / *V** Trace 
is taken over frequencies and momenta and includes fac¬ 
tor 2 accounting for the fermion spins. Transformation 
to zero temperature follows by the standard replacement 
^ _1 En ^ Jde/(2iri). 

Expanding Eq. (l42l) up to the 4-th order in (f> for T = 0 
we obtain 

S'[4>\ ~ ^sgn(/ 0 ) Y <t>q (( r o ) _1 + a 2 N$(s, x))tj>- q (43) 

Q 

V 2 V 3 094^91+92 + 93 + 94,0 i 

{ 9i } 

where the effective field self interaction is given by the 
function 

Q2i Q3i Q&) 2i — ^ ^ GpGp-\.g 1 

' V{ 1,...,4) P 

x Cp+9i+92^p-l-9i+92+93 • (44) 


Hence, we can identify the vertex of the boson-fermion 
interaction 

V B 2 = ~fom 2 B /{a 2 N) , (48) 

and the full retarded Green’s function of the boson 
D§(u, k) with the retarded self-energy 

Eg(w,fc) =/ 0 m|$(-^-, —) . (49) 

kv F PF 

Note that for very large tub that we have assumed, 
D|(w,fc) in flUD differs from D^(io,k) in (1551) only by 
a prefactor. For the attractive interaction, which we are 
now interested in (/ 0 < 0), we have m| > 0 and Eg > 0. 
For repulsive interaction instead of the bare scalar boson 
the bare vector boson would be an appropriate choice, 
cf. 371. or we may come back to the formalism given by 

Eqs m - USD- 

The spectrum follows from the solution of the Dyson 
equation 0 = [Dg(w,fc)] _1 ~ — m|(l + /o$(s, x)) for 
m-B uj,k and coincides with the solution of Eq. (HU). 
The boson spectral function is given by 


Here the sum VV runs over the 4 permutations of 4 mo¬ 
menta qi . As we will show, Eq. (l43l) is applicable also 
for /o < 0 after the replacement <fi -+ i(f> that results in 
the appearance of the prefactor sgn(/ 0 ) in Eq. (1551) . 

The general analysis of Eq. (1551) for arbitrary external 
momenta was undertaken in Ref. (.361 ]. The first term 
in Eq. (1551) (quadratic in <j >) can be interpreted as the 
inversed retarded propagator of the effective boson zero- 
sound-like mode: 


A B {ui,k) = -29tDg(w,fc), 

- 29fT p %(w, k) = V 2 A b (uj, k) . (50) 

IV. AVOIDING OF THE POMERANCHUK 
INSTABILITY BY A BOSE CONDENSATION 

A. Condensation of a scalar field for fo < —1. 


{Df) 1 W) = —sgn(/ 0 )[(To) 1 +a 2 7V$(s,x)] • (45) 

To describe modes for an attractive interaction, fo < 0, 
it is instructive to re-derive the expression for the mode 
propagator using the approach proposed in Ref. [37j. The 
local four-fermion interaction Tg can be viewed as an 
interaction induced by the exchange of a scalar heavy 
boson with the mass tob, 

Tq m B f*B,o( w i V i 
(46) 

where D B 0 is the bare retarded Green’s function of the 
heavy boson and we assumed that m B is much larger 
then typical squared frequencies and momenta, oj 2 , k 2 , in 
the problem. 


r o™4 


(w + id) 2 — k 2 — ■ 


Consider a Fermi liquid with local scalar interaction 
fo < — 1. We assume that the bosonization procedure of 
the interaction, Eqs. (1551) and (|47l) . is performed. Accord¬ 
ing to the perturbative analysis in Sect. Ill Dl for / 0 < — 1 
there are modes which grow with time. The growth of hy¬ 
drodynamic modes (first sound), cf. Eqs. (1551) and (15U1) . 
results in the formation of a mixed phase. Besides the hy¬ 
drodynamic modes the zero-sound-like modes grow with 
time, cf. Eqs. (El and EZD- We study now the oppor¬ 
tunity that the instability of the zero-sound mode may 
result in a formation of a static condensate of the scalar 
field. This leads to a decrease of the system energy and 
to a rearrangement of the excitation spectrum on the 
ground of the condensate field. 

We will exploit the simplest probing function describ¬ 
ing the complex scalar field of the form of a running wave 

ip = Woe~ iUct+i%c? , (51) 









with the condensate frequency and momentum (w c , k c ) 
and the constant amplitude tpo. The choice of the struc¬ 
ture of the order parameter is unimportant for our study 
of the stability of the system. 

Guided by the construction of the full Green’s function 
of the effective boson field, see Eq. 62D, the effective 
Lagrangian density for the condensed field (15T1) can be 
written as, cf. Ref. [22l [38], 

L = K.D3 1 (w c ,fc c )|( / ?o| 2 - \V£A{uj c ,k c )\tp 0 \ 4 . (52) 
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which corresponds to the gain in the energy density 


E b = -N 2 



e(-oj 2 (k c )), 


(59) 


where 9{x) = 1 for x > 0 and 0 otherwise. 


B. Bose condensation in a homogeneous state. 


The energy density of the condensate is given then by 
the standard relation 

Eb = oodL/duo — L . 


If the minimum of the gap Co 2 is realized at k c = 0, e.g., 
it is so for /02 > 0, the energy density Eb is minimized 
for io c = k c = 0. Then, the field amplitude and the 
condensate energy density are 


Fully equivalently, this Lagrangian density can be written 
in the form suggested by expansion (P51) . now applied for 
the running wave classical field, after the field redefinition 
m B <Po = a\J N/\f 0 \(j) 0 , 

L = \M 2 - ^A(w c , k c )\4> 0 \ 4 ■ (53) 

Here and in Eq. (Ifnil) the quantity A(w c , k c ) represents the 
self-interaction amplitude of condensed modes which cor¬ 
responds to the ring diagram with four fermion Green’s 
functions (l44l) . A(w c ,fc c ) = 6U(q c ,—q c ,q c ,~q c ). As we 
shall see below the energetically favorable is the state 
with co c = 0, and therefore we put lo c = 0 here. The 
leading order contribution to the self-interaction param¬ 
eter A(0, fc c ) as a function of the condensate momentum 
was calculated in Ref. [36| : 


4>l = N j^e{ Zf ), E B = -N 2 ^-e(z f ), (60) 

where we used the same notation Zf as in Eq. (l34l) . The 
Bose condensate is formed for f 00 < — 1. 

For the case of non-zero temperature, 7 1 / 0, in the 
mean field approximation one should just replace $ —> 
‘by, where the fermion step-function distributions are re¬ 
placed to the thermal distributions, and from Eq. CHI 
for T <C £f one recovers the critical temperature of the 
condensation T MF /eF = ^/12 z//7r valid for 0<z/Cl. 
For |/oo| » 1 from Eq. CHI one gets T MF /e f = 2|/ 0 o|/3. 

In the presence of the homogeneous condensate ( k c = 
0) setting cj) q = (j) 0 + cj)' q in (l43l) . (l53l) . varying in (f>' q and 
retaining only terms linear in <f>' q we recover the spectrum 
of overcondensate excitations 


A(0, fc c ) = 


8a 4 

1 — 

c In 

2 -f- x c 

1 — x ^ 

c In 

1 + x c 


TT 2 VpXc 

X c 

2 x c 

2x c 

1 - X c 

a 2 N 


4~x 2 J ’ 
For x c <C 1 we get 


x c = kc/pF- 


A(0, k c ) ~ a 4 A ^1 + 7 ^ 2 ^) > 


X = 


.2, ,3 


(54) 


(55) 


— l/oo|" 


$ 


/ to k 


\vpk’ pf 


A) 

VF> 


- SE V = 0, 


(61) 


where the last term arises from the interaction of excita¬ 
tions with the condensate. Making use of Eqs ED, 
and © we find 


(511^ = 2 a~ 


~ 2 N Zf + 0(w 2 , k 2 , w/k —> 0) (62) 

and derive the spectrum for Zf 1 and i<1. 


The quantity A agrees also with the result derived in 
Ref. [891 1 f° r description of the pion condensation in the 
Thomas-Fermi approximation. 

The equation of motion for the field amplitude follows 
from the variation of the action corresponding to the La¬ 
grangian density (1551) . 

—a 2 Nui 2 (k c )<fio - A(0, k c ) \4>o\ 2 <Po = 0, (56) 


Co « —i—Zfkv f . (63) 

7 r 

So, the overcondensate excitations are damped, similar to 
those we obtain for the case — 1 < /oo < 0, see Eq. (E51) . 

In the presence of the condensate (/oo < —1) the 
particle-hole interaction is described by scattering am¬ 
plitude 


where we introduce the effective boson gap 

Co 2 (k c ) = [1 - |/ o (fc c )| < i ) (0, k c )\/\fo(k c )\ 


(57) 


as it was done in the description of the pion condensation, 
see Ref. [22j. The equation for the condensate amplitude 
© has a non trivial solution for Co 2 (k c ) < 0, 


\M 2 = -n- 


Co 2 (k c 


2A ( 1 + S) 


0(-CC 2 (k c )) , (58) 


a 2 NT% = 
1 


1 

L fo 


^ ^ 1-1 

$ - v 


- /ren,0 


1 


2kvF 


a 2 N J 


+ 0(w 2 , fc 2 , to/k —> 0) 


-1 


(64) 


where in the second equality we introduce the renormal¬ 
ized local interaction 


/ren,0 ~ /ren,00 


/oo 


2/00 + 1 


(65) 
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The latter equation shows that if originally we have /oo < 
— 1 , the renormalized interaction — 1 < / ren ,oo < — 1 / 2 - 
Thus, in the presence of the boson condensate the Fermi 
liquid is free from the Pomeranchuk instability. 

Using the renormalized parameter / ren ,oo we can cal¬ 
culate how the fermion chemical potential changes in the 
presence of the condensate. We use here the standard 
expression from Ref. [ 6 j for the variation of the chemi¬ 
cal potential with the particle density in the Fermi liquid 
theory 


Sp F 


2cf,o 
3 no 


Sn + 


/ren, 00 

N 


Sn , 


( 66 ) 


here £f,o = Pp( n o)/^mp(no) . The energy density can be 
obtained from the relation E = fipdn. 

The compressibility of the system becomes now 


rr Sp F 

A = n o 

on 


+ K b 


o) 

3 mp(n 0 ) 


(! + /ren,00 ) + K b , 


(67) 


where the first term in the second line is always positive 

1 + /ren, 00 = (|/oo| ~ 1 )/( 2 |/ 00 | - 1 ) > 0 , ( 68 ) 

for /oo < — 1 under consideration. The second term is 
the condensate term: 


tv- d 2 E B 
hB=n °^ 


(69) 


Similarly, the square of the first sound velocity in the 
presence of the condensate becomes 


u 2 = 


dP 


dp 2 >m,pm 


Pf 


T (l+/ren,00) >0, (70) 


whereas in the absence of the condensate it was oc (1 + 
/oo) < 0, cf. Eq. (pH]). 

Recall that u 2 is negative within the spinodal region 
which exists, if the pressure has a van-der-Waals form. 
Thus, the formation of the scalar Bose condensate, sug¬ 
gested here, might compete with the development of the 
spinodal instability at the liquid-gas phase transition in 
Fermi liquids. From general principles, one can expect 
that the system will develop a (stable or metastable) 
condensate state at least, if a surface tension between 
liquid and gas regions is sufficiently high. Then aerosol¬ 
like mixture appearing in the course of the development 
of the spinodal decomposition will evolve more slowly 
compared to the process of the formation of the scalar 
condensate. 


C. Bose condensation in an inhomogeneous state. 

We discuss now how the momentum dependence of the 
coupling constant /o, as given by Eq. (1751) . may influ¬ 
ence the condensation of the scalar mode. The condi¬ 
tion for the appearance of a condensate with w c = 0 



FIG. 4: Phase diagram of a Fermi liquid with a scalar interac¬ 
tion in the particle hole channel described by the momentum 
dependent coupling constant /o(A:/pf) = /oo + /o 2 & 2 /p|. De¬ 
pending on the values of parameters /oo and /02 the Fermi 
liquid can be either in a stable state or in states with homo¬ 
geneous or inhomogeneous static Bose condensates (BC) of 
scalar modes. Possibility of the development of the spinodal 
instability is suppressed. 


and a finite momentum k c are determined by relations 
£b 2 (k c ) < 0, jj:Cj 2 (k c ) = 0. The critical condition for 
the appearance of the condensate can be deduced from 
the expansion of \f 0 \ui 2 (k) in (1571) for x 2 ~ x 


1 + fo(x)®( 0 ,x) 


/02 


/oo 

20 


12 


1 + /oo + 3 


(/02 - #) 2 

(/02 + ^) 

/02 + ^' 


(71) 


For /oo < — 1 this function has always a negative mini¬ 
mum at x = 0 , the second and deeper minimum appears 
at finite x, if 


/o2</o c ^ = /oo/12. (72) 

For — 1 < /oo < 0 the minimum at finite momentum can 
be realized if 


1 r 


/02 < fo 2 ~ — 7 ; 1 + 5/00 + \/l + -r/oo — 7/1 


61 


/oo 


and 


x 2 m (f 02 ) =6 


1 + /oo + \J 1 + f/oo - f/oo 
1 + ^/l - §/oo - f /oo + Ifoo 


(73) 


(74) 


The solid lines in Fig. [4] on the (/oo,/o 2 ) plane show 
the critical values /q£ given by Eqs. (1721) and (1731) . They 
form the phase diagram of the Fermi liquid with the mo¬ 
mentum dependent interaction in the scalar channel. The 
dash-dotted line shows solution of equation fo 2 Xm = /oo- 
In the region restricted by this dash-dotted line and the 
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dash continuation of the solid line we have |/oo| < |/o 2 |£m 
and the expansion (E^l) for / 0 (fc) becomes questionable. 

Near the minimum the boson gap w 2 (fc) can be pre¬ 
sented as 


u 2 (k) =u>l(k c ) +^(k 2 -k 2 ) 2 , (75) 

for 7 = const. The critical point is found from the con¬ 
dition Wj(fc c = kcr) = 0 . Beyond the critical point the 
actual value of k c follows from the minimum of the energy 
density (l59l) as a function of k c . 


V. CONDENSATE OF BOSE EXCITATIONS 
WITH NON-ZERO MOMENTUM IN A MOVING 
FERMI LIQUID 

Let us apply the constructed above formalism to the 
analysis of a possibility of condensation of zero-sound-like 
excitations with a non-zero momentum and frequency in 
a moving Fermi liquid. The main idea was formulated in 
Refs. |9l4l5|. When a medium moves in straight line with 
the velocity u > ztL = min[o;(fc)/fc] (where the minimum 
is realized at k = fco ^ 0 ), it may become energetically 
favorable to transfer a part of the momentum from the 
particles of the moving medium to a condensate of col¬ 
lective Bose excitations with the momentum kg. The 
condensation may occur, if in the spectrum branch oj(k) 
there is a region with a small energy at sufficiently large 
momenta. 

As in Ref. 0 , we consider a fluid element of the 
medium with the mass density p moving with a non- 
relativistic constant velocity u. The quasiparticle energy 
w(fc) in the rest frame of the fluid is determined from the 
dispersion relation 


^D~\uj,k) = 0. (76) 

We continue to exploit the complex scalar condensate 
field described by the simplest running-wave probing 
function, cf. Eq. m, and the Lagrangian density (1551) . 
but now for the condensate of excitations. 

The appearance of the condensate with a finite mo¬ 
mentum fco, frequency w = w(fco) and an amplitude tpo 
leads to a change of the fluid velocity from u to Mfi n , as 
it is required by the momentum conservation 

pu = pufi n + fc 0 ^o" 1 |<( , ol > ( 77 ) 

where A?o^o~ 1 | < ^ ) ol the density of the momentum of the 
condensate of the boson quasiparticles with the quasipar¬ 
ticle weight 


z o\ko) = [^StD^fakj 


u(ko),k 0 


> 0 . 


(78) 


which takes a part of the momentum, the energy density 
becomes 


Efin = \pu\ n +U](k 0 )z 0 1 \(j) o | 2 + ±A(w(fc 0 ),fc 0 )|^o| 4 - 

(79) 


Here the last two terms appear because of the classical 
field of the condensate of excitations. The gain in the 
energy density due to the condensation, SE = Ea n — E- m , 
is equal to 

5E = -[ufc,o - w(fc 0 )] ZQ 1 (k 0 ) \cj) 0 \ 2 + iA|0 o | 4 , (80) 

where 


A = A(w(fc 0 ), k 0 ) + {Z^\k 0 )) 2 k 2 Jp . (81) 

For oj = 0, A(0,fco) is calculated explicitly, cf. Eq. (l55l) . 
Note that above equations hold also for A = 0. 

The condensate of excitations is generated for the ve¬ 
locity of the medium exceeding the Landau critical ve¬ 
locity, u > ml = w(fco)/fc’o, where the direction of the 
condensate vector fc 0 coincides with the direction of the 
velocity, fc 0 || u, and the magnitude fco is determined by 
the equation w(fco)/fco = doj(fco)/dfc. The gain in the 
energy density after the formation of the classical con¬ 
densate field with the amplitude <j >o and the momentum 
fco is then 

SE = -ZQ 1 (k 0 )[uk 0 - w(fco)]^o + ^A^o ■ (82) 

The amplitude of the condensate field is found by mini¬ 
mization of the energy. From (18211 one gets 

<Po — z o (fc'o) ~ S(U — lie). 

The resulting velocity of the medium becomes 
(u - u c )9{u - u c ) 

hn c l + [Z 0 - 1 (fc 0 )] 2 fc 0 2 /(A( w (fco),fco)p 

For a small A, we have ua n = u c + O(A). 

For the repulsive interaction /o > 0 there is real zero- 
sound branch of excitations w s (fc) ~ ki’F(so + S 2 x 2 + 
S 4 :r 4 ), where the parameters s, depend on the cou¬ 
pling constants /oo and /02 according to Eqs. m, <m, 
and m- As shown in Sect. nm the ratio ui s [k)/(kv f) 
has a minimum at fco = Pf\/— S2/(2 s4 ) provided /02 is 
smaller than / crit ,02 = —/oo/[ 12 ( s o _ l) 2 ]- The Lan¬ 
dau critical velocity of the medium is equal to ml/mf = 
Wg(fco) / (^Ffco) ~ so — s|/(2s 4 ). The quasiparticle weight 
of the zero-sound mode CHI) is 


(83) 

(84) 


z o\k 0 ) 


a 2 N d$(s,x) 
kgVF 0s 


^s(fcp) fcp 
VpkQ ’ Pp 


(85) 


The amplitude of the condensate field (l 88 l) can be written 
as 


Pq — Zq (fco)fco- 


-9{u - m l ) ■ 


If in the absence of the condensate of excitations the 
energy density of the liquid element was E- ln = pu 2 /2 , 
then in the presence of the condensate of excitations, 


A 


( 86 ) 
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The energy density gained owing to the condensation of 
the excitations is 


SE = -k 2 0 (Z~ 1 (k 0 )) 


2 (u-u h ) 2 


2A 


9{u — ul). (87) 


For a small A, Eqs. HMD, dm) simplify as 


2 (u-u h ) 

SE =-^(u - u L ) 2 9(u - u L ). (88) 


VI. CONCLUDING REMARKS 

We described the spectrum of scalar excitations in nor¬ 
mal Fermi liquids for various values of the Landau param¬ 
eter /o in the particle-hole channel and for different mod¬ 
els of its momentum dependence. For /o > 0 we found a 
condition on the momentum dependence of fo(k) when 
the zero-sound excitations with a non-zero momentum 
can be produced in the medium moving with the veloc¬ 
ity larger than the Landau critical velocity. Such ex¬ 
citations will form an inhomogeneous Bose condensate. 
For — 1 < f 0 < 0 there exist only damped excitations. 
For / 0 < — 1 we studied the instability of the spectrum 
with respect to the growth of the zero-sound-like exci¬ 
tations (Pomeranchuck instability) and the excitations 
of the first sound. The surface tension coefficient above 
which the zero-sound-like mode grows more rapidly than 
the hydrodynamic one (for ideal hydrodynamics) is es¬ 
timated. Then we derived an effective Lagrangian for 


the zero-sound-like modes by performing bosonization of 
a local fermion-fermion interaction. We argue that the 
Fermi liquid with fo < — 1 might become stable owing 
to appearance of the static (homogeneous or inhomoge¬ 
neous) Bose condensate of the scalar quanta. Properties 
of the novel condensate state are investigated. 

In the future it would be important to search for a pos¬ 
sibility of realization of this phenomenon in some Fermi 
liquids. It would be interesting to perform a careful com¬ 
parative study of the possibilities of the Bose condensa¬ 
tion and the ordinary spinodal instability appearing at 
the first-order phase transitions in the Fermi systems. 
For this, an explicit expression for the energy density 
functional for a specific system, e.g. for the nuclear mat¬ 
ter, is required. Such a programme might be realized 
within the relativistic mean-field model. Similarly, we 
expect a stabilization of the Fermi liquid with a strong 
attractive spin-spin interaction go < —1 by a conden¬ 
sate of a virtual boson field at certain conditions. These 
questions will be considered elsewhere. 
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